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We theoretically investigate time-dependent resonant tunneling via two discrete states in an ex- 
perimentally relevant setup. Our results show that the dc transport through the system can be 
controled by applying external irradiation with a frequency which matches the energy difference 
between the discrete states. We predict resonant phenomena which should be easily observable in 
experiments. 



Time-dependent tunneling phenomena have received 
increasing attention in recent years. In the early eighties, 
Biittikcr and Landauer studied the tunneling tunc needed 
for an electron to traverse a potential barrier .El More re- 
cent theoretical work focused on the time-dependence 
of resonant tunneling using an effective Schrodinger 
equations and on a description of the time-dependent 
current through mesoscopic structures in terms of non- 
equilibrium Green's functions.u In addition, the consid- 
erable improvement in nanofabrication techniques facil- 
itated some interesting experimental studies. Kouwen- 
hoven et al. measured the photon-assisted tunneling cur- 
rent through a single quantum dot with an effectively 
continuous level spectrum due to thermal smearing.!] Van 
der Vaart et al. studied the dc current through a double 
dot system with well developed 0D states in each dot and 
clearly resolved resonances between energy levels in both 
dots. El The sharp resonance features make it very tempt- 
ing to perform experiments with time-dependent fields. 
The dc current through such a structure in the presence 
of oscillating fields may be expected to display interesting 
new phenomena, not observable in a single dot. 

In this paper we use of the density matrix approach of 
Ref. H in which the resonant states are described by a 
time-dependent tunneling Hamiltonian. Transitions be- 
tween non-resonant states of the system are taken into 
account through a master equation for the density matrix 
elements. We calculate the photoresponse of the system 
in several experimentally relevant limits and derive an 
explicit expression for the shape of the resonant peaks in 
the case of an external perturbation with arbirary am- 
plitude. Close to resonance the dc current is found to be 
very sensitive to the oscillating field. The satellite reso- 
nances induced by the external oscillating field can be of 
the same order of magnitude as the main static resonance 
with an even smaller width. 

The system under consideration (Fig. |l|) consists of 
two quantum dots A and B in series. The dots are con- 
nected by tunnel junctions to two large reservoirs L and 
R, which are assumed to have continuous energy level 
spectra and are filled up to their respective Fermi ener- 
gies. If we neglect all tunneling processes, a system of 
discrete many-body states is formed in each dot. The 
best conditions for transport occur when it costs no en- 
ergy to transfer an electron between the dots, i.e. the 



energy difference between a state with one extra electron 
in the left dot and a state with one extra electron in the 
right dot is zero. In the experiment in Ref. || this energy 
difference could be tuned by an external gate voltage. 
The current through the system vs. gate voltage consists 
of a series of peaks corresponding to the resonances be- 
tween different discrete states. There could be a variety 
of different transport processes occuring in a resonance 
point as described in Ref. || We concentrate on the sim- 
plest experimentally relevant case, namely when the res- 
onance occurs between the ground states of both dots. 
We assume that the bias voltage is much larger than the 
temperature and the energy difference between the states 
in resonance. Consequently, electrons can only enter the 
two-dot system from the left and leave it only to the right. 
Transitions from the left and to the right lead are possible 
with rates Tl and r#, respectively. Here and throughout 
the paper units are used such that h = 1. We will as- 
sume that it is impossible, due to the Coulomb blockade, 
for an electron to tunnel into the system while another 
electron is still present in either one of the dots. We con- 
centrate on two states only and disregard other states, 
which is allowed in the neighborhood of the resonance. 
The system can also be in a third state |0), when there is 
no extra electron in either one of the dots. The energies 
of the resonant states, which both lie well between the 
electrochemical potentials in the left and right lead, ini- 
tially differ by an amount cq. Under these conditions the 
transport through the system depends only very weakly 
on the bias voltage, but does depend strongly on the gate 
electrode via the energy difference eo. We assume that 
a time-dependent oscillating signal is applied to the gate 
electrode, so that the time-dependent energy difference 
becomes: e(t) = eo + ecostut, where e is the amplitude 
and uj the frequency of the externally applied signal. 

The dynamics of the resonant states, \a) and \b), is 
governed by the time dependent tunneling Hamiltonian 
H(t) — Ho(t) + Ht, where Ho is given by 

H (t) = le(t)(\a){a\-\b){b\) (1) 

and Ht describes the coupling between the dots that 
introduces mixing between the eigenstates \a) and \b) of 
the system: 

H T = T(\a)(b\ + \b)(a\). (2) 
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The average current through the system is given by 



The first order correction to the stationary solution is 



(I)/e=Tr(pl), 
where X is the current operator: 

I = *T(|a)(6|-|6)(a|) 



(3) 



(4) 



and p is the density matrix for the two-level system. We 
describe transitions between different states in the den- 
sity matrix approach.tl The equations for the density ma- 
trix elements read 



dp a 
dt 
dp b 
dt 
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dt 

dpba 

dt 
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r« Pab + te(t) Pab + iT {p b - Pa) 
^ Pba - ie(t) Pba - iT {pb ~ Pa), 



(5a) 
(5b) 
(5c) 
(5d) 



where p a , pb and p a = 1 — p a — pb denote the probabili- 
ties for an electron to be in the left dot, the right dot, or 
in neither dot, respectively, and p a b — p% a are the non- 
diagonal density matrix elements. In these equations the 
terms proportional to and Tr describe the transitions 
to and from the reservoirs between the states |0) and \a) 
and the states \b) and |0), respectively. All other terms 
follow from the Liouville equation: idp/dt = [H, p]. Note 
that the rates Tl and do not enter the equations in 
a symmetric way. Tr describes the decay of the reso- 
nant states whereas Tl describes the build-up of these 
resonant states. 

The relevant energy scales of the system are the tran- 
sition rates Tl and Tr, the tunneling amplitude T, and 
the frequency uj and amplitude e of the applied pertur- 
bation. There are three limiting cases for which we can 
develop an analytical approach to the problem. They are 
complementary and essentially cover all the interesting 
physics. 

We will first consider the limiting case of a small per- 
turbation amplitude; e -C u>, T,Tl^r. Using the fact that 
po = 1 — p a — pb, we rewrite Eqs. (p[) in matrix notation: 



dp 
dt 



f + f + 6q + £ COS Ujt) p + C, 



(6) 



where p = (pa, Pb, Pab, Pba) T , c = (r L ,0,0,0) T and the 
matrices f, T, £q and e correspond to Eqs. (|^). The 
stationary solution of these equations without irradiation 



is 



Po 



-(T + T + eo) 



(7) 



This determines the shape of the stationary resonant 
peaks observed by Van der Vaart et al.-Q 



-^stat 



T2(2 + r fl /r L ) + r|/4 + e2 
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(8) 



pi = p^ cxv(iujt) + p x cxp(-icut) , 

with the positive and negative frequency part 
spectively: 



(9) 



rc- 



Pl = -(T + T+ e T wl)' 1 -^, 



(10) 



/ being the unit matrix. This contribution contains only 
oscillatory terms, which average out when calculating the 
dc current. We therefore also determine the second order 
correction terms (proportional to e 2 ) and obtain 



Pi 



-(r + T + 6 ) 



-\{Pi + + Pi)- 



(ii) 



Using the density matrix elements pi we may calculate 
the photoresponse of the system. This quantity can be 
easily measured experimentally by slowly modulating the 
irradiation amplitude.cl In Fig. a plot is given of the 
photoresponse as a function of eo and u> for = Tr ■ — 
0.2T. The figure clearly shows resonant satellite peaks for 
uj and eo satisfying uj 2 = €q + AT 2 , i.e. resonant modes 
occur when the frequency of the applied per turbation 
matches the renormalized energy difference yel-l- 4T 2 
of the two levels. For frequencies below 2T there are no 
satellite peaks because the energy Ticu of the photon is 
smaller than the energy level spacing. The evolution of a 
resonant satellite peak is shown in Fig. [?], where a current 
peak for u> = iT has been plotted vs. eo for different 
values of Tr/T = Tl/T. We see that the peak can be 
seen even at moderately large values of Tr/T, but the 
best resonance conditions occur when Y r <C max(T, eo) 
and u = ^e 2 +4T 2 . 

We have developed a second approach which allows us 
to explore the satellite peak at arbitrary values of irradi- 
ation amplitude e under the conditions mentioned above. 
Substituting p = po + p+(t) exp(icjf) + p~(t) exp(—iojt) in 
Eq. ([|) and neglecting terms proportional to exp(±2iw£) 
we obtain: 

^ - f pb + (eo + f)p Q + | (p+ +p-) + c, (12a) 

^=fp + + (e Q + f-iuji)p + + ^po> (12b) 

^ = fp_ + (e + f + iui)p- + ~p . (12c) 

Near the resonance point we can approximate the so- 
lution p by an expansion in terms of the eigenvectors 
of the matrix eo + T: po — a\V\ + U2V2, P+ = a +v+ 
and p_ = where v\ and V2 are the eigenvec- 

tors with eigenvalue and v± those with eigenvalues 
±iyeo + 4T 2 . We obtain a set of four closed equations 
for the coef ficients 0^2,+,- by taki ng th e inner product 
of E q. ( | 12a| ) with V\ and V2, of Eq. (12b) with v+, and of 



Eq. (12c ) with w_ . Solving for the stationary solution and 
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calculating the current profile near the resonance point 
e r = y/uj 2 — 4T 2 results in a Lorentzian line shape: 



with height 



I/e 



In 



W 2 + (e - e r ) 2 ' 
~e 2 Y R (a 2 - 4) 



and half- width at half maximum 



2y/a 2 - 4 
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(13) 



(14) 



(15) 



where a = u/T, (3 = lVr L + 2, and 7 = a 2 + (3 - 4. 

In the limit of small amplitude the height of the cur- 
rent peak scales with the square of e: /max = e 2 (a: 2 — 
4)/7 2 r^, whereas the width remains constant: w = 
iar^/i/a 2 — 4 ~ Fr, consistent with the results pre- 
sented in Fig. [| At a — 2 (corresponding to eo = 0) the 
peak vanishes, as seen in Fig. g. With further increase of 
e the current saturates at a value of I sa t = T n(a 2 — 4) / P~y 
which is of the order of the height of the stationary peak. 
This saturation occurs at relatively small e ~ Tr. The 
width of the peak increases with growing e. At e Tr 
it is proportional to e. 

Therefore we have shown that, under good resonance 
conditions, the current is very sensitive to the external 
irradiation. A relatively weak irradiation induces a big 
satellite peak that has a much smaller width than the 
stationary one. 

For small tunneling amplitudes T, provided eo 3> I\r, 
the height scales with T 2 : I max = T 2 € 2 /Truj 2 , and the 
half-width reduces tow = jFr. These results agree with 
the expression for the photon-assisted tunneling current 
derived below, where we consider our third approach, in 
which the tunneling amplitude is small compared to all 
other energy scales in the system; T< e,w, r^^. 

First we perform a transformation on the density ma- 
trix which leaves the diagonal elements invariant and 
which changes the non-diagonal elements as follows: 



Pab = Pab exp \-i 



dr e(r) 



(16) 



This transformation eliminates the explicit time depen- 
dence in Eqs. (0c ) and (||d) and introduces it in the trans- 
formed tunneling amplitude. The equations for the non- 
diagonal density matrix elements now assume the form 



dpab 

dt 



te 



with time-dependent tunneling amplitude 



T(t) = Texp [i / dr e(r) 



T R ) p ab + if(t){p b - p a ), (17) 



(18) 



The equations for p\, a are simply the complex conjugate 
of Eqs. ([l6l), ( p7| ) and (|is|). For p ab ^ in the lowest non- 
vanishing order in T we obtain: 



2 r _R + i{ui - eo) 
Expanding T(t) in a Fourier series; 



(19) 



T{t)=T Me/u)e*p(-inuJt), (20) 

n— — 00 

and calculating the dc current, Eq. (||), results in: 

(I)/e = T 2 T R £ lr2 f, (£>) r- 2 , (21) 

where J n is a n'th order Bessel function of the first kind. 
This equation for the current is similar to the expression 
found by Tien and GordonEI for the photon assisted tun- 
neling current through a superconducting tunnel junc- 
tion. Note, however, that in the Tien-Gordon case the 
current has been considered as a function of bias voltage 
whereas in our case it is a function of eo , the energy shift 
of the discrete levels. Analogously, the alternating field 
is not applied in the bias direction but rather to the gate 
electrodes. 

In Fig. |] the current has been plotted as a function 
of eo and e. The figure clearly shows that the current 
is composed of a number of satellite peaks each sepa- 
rated by the photon energy Tim. With increasing am- 
plitude e, the number of visible current peaks increases. 
The peaks all have the same width Tr and have heights 
given by 4T 2 J 2 (£/u})/Tr. In the limit of small ampli- 
tude the height of the n = 1 satellite peak reduces to 
Imax = T 2 e 2 /Truj 2 , identical to our earlier result. Note 
that Eq. (^lj) for the current no longer contains T^. Be- 
cause the tunnel rate from lead to dot is much larger than 
the tunnel rate between the dots, the width of the level 
is in this case determined by Tr and T only. 

In conclusion, we have presented a complete theoret- 
ical picture of the dc transport through a double quan- 
tum dot in the presence of external harmonic irradia- 
tion. The photoresponse of the system exhibits extra 
resonant peaks when the frequency of the external irradi- 
ation matches the energy difference between the discrete 
states. At further increase of the irradiation intensity 
this satellite peak becomes of the same order of magni- 
tude as the main peak, but preserves the much smaller 
width. At small tunneling amplitudes and large irradia- 
tion amplitude extra satellite peaks appear in a pattern 
similar to that obtained for a Josephson junction by Tien 
and Gordon.Q 
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FIG. 1. Schematic picture of the system. Two quantum 
dots A and B are coupled to leads L and R via tunnel junc- 
tions. Transitions are possible with rates Fl and T_r. The 
tunneling rate between the dots is T and the energy differ- 
ence between the levels is denoted by e(t). 

FIG. 2. Scaled photoresponse of the system as a function 
of the energy difference eo/T between the levels and the fre- 
quency w/T of the applied signal. The plot was made with 
T L = r R = 0.2T. 

FIG. 3. Evolution of a satellite peak for ratios Tl/T= 
Tr/T= 0.2, 0.5, 1.0 and 3.0 respectively. The plots were 
made for a frequency of lj — 3T. 



FIG. 4. Scaled current through the dots as a function of the 
energy difference eo/T between the levels and the amplitude 
e/T of the applied signal. The plot was made with Tr = T 
and w = 10T. 
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